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After a brief summary of the neutrino oscillation formalism and the solar 
neutrino sources and experiments I discuss the matter effect on solar neutrino 
oscillation. Then I discuss how the resulting alternative solutions are exper- 
imentally resolved in favour of the LMA solution, with particular exphasis 
on the SK, SNO and KL data. 
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The last four years have been widely described as the golden years of solar 
neutrino physics, thanks to three pioneering experiments - SuperKamiokande 
(SK), Solar Neutrino Observatory (SNO) and KaniLAND (KL). They have 
provided for the first time a unique solution to the solar neutrino anomaly 
in terms of neutrino oscillation with unambiguous mass and mixing parame- 
ters. I shall give an overview of the subject with particular emphasis on the 
role of these experiments. After a brief summary of the neutrino oscillation 
formalism and the solar neutrino sources and experiments 1 shall discuss the 
matter effect on solar neutrino oscillation. We shall see how it leads to four 
alternative solutions to the solar neutrino anomaly and then their experimen- 
tal resolution in favour of the so called Large Mixing Angle (LMA) solution 
over the last four years. 

Neutrino Mixing and Oscillation: It was already noted by Pontecorvo 
back in the sixties that if the neutrinos have non-zero mass then there will 

in general be mixing between the flavour and the mass eigenstates, which 
will lead to neutrino oscillation [1]. We assume for simplicity two neutrino 
flavours, in which case mixing can be discribed by one angle 9, i.e. 

In fact it provides a very good approximation to the three neutrino mixing 
scenario for solar neutrino oscillation, with the second neutrino representing 
a mixture of the i/^ and fr flavours. Now each mass eigenstate propagates 
with its own phase 

^-iiEt-pe)^^-i^^ (2) 

where we have made the relativistic approximation, E — p + m^/2p, since 
neutrino masses are much smaller than their kinetic energy. Thus a pro- 
duced at the origin will propagate as 

2 2 

Vf. vi cos9e~^ + 1^2 sin 9e^^ . (3) 

Decomposing the i'i^2 back into fe,^ after a distance £ one sees that the z/^ 
terms do not cancel, which implies neutrino oscillation. In fact the coefficient 
of the term represents the probability of — > oscillation, i.e. 
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= sin2 2^sm^(^£), (4) 

where the first factor represents the amplitude and the second factor the 
phase of oscillation, with Am^ = ml — m\. Converting to convenient units 
for Am^(eV'^), £(m) and E (MeV) one gets 

p^^^^^(£) = sin2 2^sin2(1.3Am2£/E), (5) 

which corresponds to an oscillation wavelength 

^ = ■ ^ 2.4E/Am2. (6) 
1.3 Am^ 

Thus for large mixing angle (sin^ 20 ~ 1) one expects the following pattern 
of oscillation probability from eqs. (0) and ©, where the factor of \ in the 
last case comes from averaging over the phase factor. 

I < A = A/2 > A 

(7) 

Pu^^u, sin2 20-1 \ sin^ 20 ~ | 

Note that the corresponding survival probability is given by the remainder, 
i.e. 

Pee = Pve^Ui: ^ Pye—*^^^- (S) 

Now the typical energy of coming from a nuclear reactor or the sun is ~ 1 
MeV. The typical distance between the source and the dector is a few hundred 
Km {I ~ 10^ m) for the long base line KamLAND reactor experiment, while 
^ ~ 10^^ m for solar neutrino experiments. Thus one sees from eqs. © 
and ^ that the KamLAND and the solar neutrino experiments can probe 
neutrino mass down to Am^ ~ 10~^ eV^ and 10~^^ eV^ respectively, which 
is far beyond the reach of any other method of mass measurement. 

Solar Neutrino Sources and Experiments: The main sources of solar 
neutrinos are the pp chains of nuclear reaction, taking place at the solar core, 
which convert protons into ^He {a particle). They are 

(I) pp ^ + e+ + i^e, "^H + p-^ ^He + 7, ^He + ^He ^He + 2p; or 

(II) ^He + ^He ^ Be + 7, ^Be + e~ ^ ^Li + z/^, '^Li+p^2 ^He; or 
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(Ill) ^Be+p^^B + 7, ^B ^Be* + e+ + Ue, ^Be* 2 "^He. 

While most of this conversion takes place by the shortest path (I) a small 
fraction (15%) follows a detour via '^Be (II); and a tiny fraction (0.1%) of 
the latter follows a still longer detour via ^B (III). The resulting neutrinos 
are (1) the low energy pp neutrino, (II) the intermediate energy Be neutrino 
and (III) the relatively high energy B neutrino, with decreasing order of flux. 
The standard solar model (SSM) prediction for these fluxes is shown in Fig. 
1 from BP 2000 [2] . It also shows the neutrino energy ranges covered by the 
different solar neutrino experiments. 

The Gallium [3] and the Chlorine [4] experiments are based on the charged 
current reactions 

Ue + ^^Ga^e- + ^^Ge and i/g + ^^CZ ^ + =^Mr. (9) 

The produced ^^Ge and ^'^Ar are periodically extracted by radiochemical 
method, from which the incident neutrino fluxes are estimated. The value of 
measured flux R relative to the SSM prediction gives the survival proba- 
bility Pee- The SK is a real-time water Cerenkov experiment [5], based on the 

elastic scattering of neutrino on electron. The elastic scattering is dominated 

cc 

by the charged current reaction z/g + e~ — > u^ + e^ , but it also has a limited 
sensitivity to the neutral current reaction i/g,/* + e~ — > t'e,/* + e~. Thus 

a^^ 1 

= ^ee + ^r^c+cd ^ " " + g(l " Pee)- (10) 

This experiment can also measure the energy and direction of the incident 
neutrino from those of the outgoing electron. The SNO is a Cerenkov ex- 
periment with a heavy water target, which can detect both the charged and 
neutral current events [6-9] 

Ue + d p + p + e~, (11) 

i^e,/. + d ^ Ue,^+p + n. (12) 

In the 1st phase of the experiment the NC events were detected via neutron 
capture on dcutcron, n + c? ^ t + 7 [7]. In the 2nd phase salt was added to 
the heavy water target to enhance the NC detection efficiency via neutron 
capture on chlorine, n + ^^Cl — > ^^Cl -|- 7 [8]. In the 3rd phase of this 
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experiment going on now ^He gas filled counters are inserted into the heavy 
water target to detect NC events via n + ^He — > t + p [9]. 

Table 1 shows the energy thresholds of the above four experiments along 
with the compositions of the corresponding solar neutrino spectra. It also 
shows the corresponding survival probability Pee measured by the rates of 
the charged current reactions relative to the SSM prediction [2]. For the SK 
experiment the survival probability calculated from i?e^ via eq. (jlUj) is shown 
in parentheses. It shows that the survival probability is slightly above 1/2 
for low energy z/g, falling to 1/3 at high energy. To understand its magnitude 
and the energy dependence we have to consider the effect of solar matter on 
neutrino oscillation. 

Table 1. The z/g survival probability Pee measured by the CC event rates of 
various solar neutrino experiments relative to the SSM prediction. For SK the 
Pee obtained after NC correction is shown in bracket. The energy threshold 
and composition of neutrino beam are also shown for each experiment. 



Experiment 


Gallium 


Chlorine 


SK 


SNO-I 


R 


0.55 ±0.03 


0.33 ±0.03 


0.465 ± 0.015 
(0.36 ±0.015) 


0.35 ± 0.03 


Eth (MeV) 


0.2 


0.8 


5 


5 


Composition 


pp (55%), 
Be (25%), B (10%) 


B (75%), 
Be (15%) 


B (100%) 


B (100%) 



Matter Enhancement (Resonant Conversion): Propagation through 
solar matter gives an induced mass to z/g, which can have profound effect on 
neutrino oscillation. This is known as MSW effect after its authors [10]. It 
arises from the charged current interaction of z/e with solar electrons, while the 
neutral current interaction has no net effect since it is flavour independent. 
Adding this interaction energy density to the free particle wave equation 
gives 

where Hint = V^GNe{0) for with G and denoting Fermi cou- 

pling and solar electron density. The quantity on the rhs numerator can be 
regarded as an effective mass or energy, i.e. 

-'^^ u :)("o' ,::i)(r;)^c^r= 
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2y/2EGNe -scAm^ 
—scAm? c^m\ + s^mj 



(14) 



where s, c denote sin 6, cos 6'. 

To get a simple picture let us assume for the moment that sin^ -C 1, 
so that the nondiagonal elements are small. Then we can identify the two 
diagonal elements with the eigenvalues and the corresponding eigenstates 
z/i^2 with the flavour eigenstates z/e,^ respectively. Fig. 2 shows the two 
eigenvalues Ai,2 against the solar electron density. At the solar surface (iVe = 
0) the 1st eigenvalue (m^) is smaller than the 2nd But Ai increases 

steadily with and becomes much larger than A2 at the solar core. The 
cross-over occurs at a critical density 

corresponding to = Mg^. Note however that the two eigenvalues actually 
never cross. There is a minimum gap between them given by the nondiagonal 
element, i.e. 

r = Am^ sin^ 26. (16) 

This means that a produced at the solar core will come out as V2, provided 
the transition probability between the two energy levels remains small. 

It is easy to show that this remarkable result does not depend on the 
sin 6* <^ 1 assumption. Consider the effective mixing angle 6m in matter, 
which diagonalises the above matrix, i.e. 

2Mg _ sin 2^ 

tan2^M - 1^,2 _ ^,2| - I ^^^2^ _ 2^GENjAm^\ ' ^ 

The electron density at the solar core is ^ N^, so that the 2nd term 
in the denominator of eq. ()17|) is much larger than the 1st. This means 
6m ^ 1 at the solar core for any vacuum mixing angle 6\ so that the u^. 
produced there is dominated by the vi component. At the critical density 
the denominator of eq. (fTTj) vanishes, which corresponds to maximal mixing 
between the two components, again for any value of 6. This is why it is called 
matter enhanced (or resonant) conversion. It comes out from the sun as V2 
with 

Pee = sin^ 6, (18) 
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provided the transition probability between the two levels remains small 
through out the propagation and in particular in the critical density region. 
This transition probability is given by the Landau-Zenner formula, i.e. 

\{d\,/M\ dNjdi 
T = e ^f, ^ = oc A^^^ — ^, (19) 

where A represents the oscillation wavelength in matter (which should not 
be confused with the eigenvalues Ai,2)- If the solar electron density varies so 
slowly that the resulting variation in the 1st eigenvalue over an oscillation 
wavelength is small compared to the gap between the two, then 7^1 and 
the transition rate is exponentially suppressed. This is called the adiabatic 
condition. Thus the two conditions for the solar to emerge as z/2 can be 
written as 

Am^ cos 261 ^ Am^ sin^ 26' 

2V2GN^ ^ ^ 2 cos 29 {dNe/dme)c' ^ ^ 

where the 1st inequality ensures > andthe 2nd one is the adiabatic 
condition. 

Fig. 3 shows the triangular region in the Am^ — sin^ 29 plot satisfying 
the above two conditions for a typical neutrino energy of ii^ = 1 MeV. The 
horizontal side follows from the 1st inequality, which gives a practically con- 
stant upper limit of Am^ since cos 2^ ~ 1. The 2nd inequality (adiabatic 
condition) gives a lower limit on sin^ 29. Moreover since this condition im- 
plies a lower limit on the product Am^ sin^ 29, it corresponds to the diagonal 
line on the log- log plot. The vertical side is simply given by the physical 
boundary, corresponding to maximal mixing. This is the socalled MSW tri- 
angle. The indicated survival probabilities out side the triangle follow from 
the vacuum oscillation formulae (|7j) and (jH)), while that inside is given by 
eq. (HHI). Thus P^e < 1/2 inside the MSW triangle and > 1/2 outside it, 
except for the oscillation maximum at the bottom (£ ~ A/2) where Pee goes 
down to cos^ 29. Finally the earth matter effect gives a small but positive 
regeneration probability, which implies a day-night a symmetry - i.e. the sun 
shines a little brighter at night in the z/g beam. After day-night averaging 
one expects a z/g regeneration probability 

. VEsm^29 _ 0.66 / AmV^ \ 

^'^^ 4(l-2r/£;cos2^ + r7i)' pY, [lO'^-^ eV ^ ' ^ ^ 

where p is matter density in the earth in gm / cc and Yg the average number of 
electrons per nucleon. For favourable values of Am^ and 9 6Preg can go upto 
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0.15 as indicated in the Figure. It is important to note that the positions of 
the MSW triangle, the earth regeneration region and the vacuum oscillation 
maximum depend only on the ratio Am^/i?, as one can see from the relevant 
formulae. Thus their positions on the right hand scale hold at all energies. 

Four Alternative Solutions: Fig. 3 marks four regions in the mass 
and mixing parameter space, which can explain the magnitude and energy 
dependence of the survival probability Pee shown in Table 1. They cor- 
respond to the socalled Large Mixing Angle (LMA), Small Mixing Angle 
(SMA), Low Mass (LOW) and Vacuum Oscillation (VAC) solutions. For the 
LMA and SMA solutions (Am^ ~ 10"^ eV^) the low energy Ga experiment 
(£" < 1 MeV) falls above the MSW triangle in Am'^/E, while the SK and 
SNO experiments {E ^ 1 MeV) fall inside it. Therefore the solar matter 
effect can explain the observed decrease of the survial probability with in- 
creasing energy. For the LOW solution the low energy Ga experiment is 
pushed up to the region indicated by the dashed line, where it gets an ad- 
ditional contribution to the survival probability from the earth regeneration 
effect. Finally the VAC solution explains the energy dependence of the sur- 
vival probability via the energy dependence of the oscillation phase in eq. 
(jSj). Fig. 4 shows the predicted survival probabilities of the four solutions as 
functions of neutrino energy. The LMA and LOW solutions predict mild and 
monotonic energy dependence, while the SMA and VAC solutions predict 
very strong and nonmonotonic energy dependence. 

Experimental Resolution in Favour of LMA: The survival rates in 
Table 1 show a slight preference for a nonmonotonic energy dependence, since 
the intermediate energy Chlorine experiment shows a little lower rate than 
SK and SNO. Therefore the SMA and VAC were the favoured solutions in the 
early days. However the situation changed drastically with the measurement 
of the energy spectrum by SK [5], shown in Fig. 5. It shows very little 
energy depence in clear disagreement with the predictions of the SMA and 
VAC solutions in Fig. 4. In particular the SMA solution in Fig. 4 can 
not reconcile a low survival rate of Pee — 1/3 with an energy independent 
spectrum, in conflict with the SK data. The first charged current data from 
SNO [6] agreed with the low survival rate as well as the energy independent 
spectrum of the SK data. Thus the global analyses of the solar neutrino data 
at this stage ruled out the SMA and VAC solutions in favour of LMA and 
LOW [11]. 
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Then came the first neutral current data from SNO [7]. Being flavour 
independent the neutral current reaction is unaffected by neutrino oscillation. 
Hence it can be used to measure the Boron neutrino flux, for which the SSM 
has a large uncertainty (see Fig. 1). The SNO neutral current measurement 
of this flux was in agreement with the SSM prediction and signiflcantly more 
precise than the latter. Fig. 6 shows the results of global flt with and without 
the SNO (NC) data [12]. The left and the middle panels show two different 
methods of using the NC data and give very similar results. They strongly 
favour the LMA solution while barely allowing LOW at 3a level. The reason 
is that the earth regeneration contribution is too small to account for the 
large survival rate of the Ga compared to those of the SK and SNO (CC) 
experiments. Without the SNO (NC) data it was possible to effectively push 
up the survival rate of the latter experiments from 0.35 to 0.45 by exploiting 
the large uncertainty in the Boron neutrino flux, and hence accommodate the 
LOW solution as shown in the right pannel (see also Fig. 4). However this 
was no longer possible after the SNO (NC) data [12,13] due to the improved 
precission of the Boron neutrino flux. 

Confirmation and Shcirpening of the LMA Solution: Independent 
conflrmation of the LMA solution came from the reactor antineutrino data 
of the KamLAND experiment [14]. assuming CPT invariance. It is a 1 kiloton 
liquid scintillator experiment detecting z/g from the Japanese nuclear reactors 
via 

z7g + p ^ e+ + n. (22) 

It also measures the incident i7e energy E via the visible scintillation energy 
produced by the positron and its annihilation with a target electron, i.e. 

E^is^ E + me + mp-mn^ E - 0.8 MeV. (23) 

The mean base hne distance of the detector from the reactors is {£) ~ 180 km, 
which means it is sensitive to the Am^^lO"^ eV^ region as mentioned earlier. 
Thus the experiment was designed to probe the LMA region. It was shown 
in [15] that if the survival rate seen at KL is < 0.9, then it would rule out 
the LOW solution at 3a level. The flrst KL result from 162ty data showed a 
survival rate R — 0.611 ±.085 ±.041 [14]. This was in perfect agreement with 
the LMA prediction and ruled out LOW at 5a level. Moreover the observed 
spectral distortion of the KL data, taken together with global solar neutrino 
data, conflned the LMA solution to two subregions around Am^ — 7 and 
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14 X 10~^ eV^, as shown in Fig. 7 [16]. They correspond to the 1st and 2nd 
oscillation minima ((^) = A and 2A) for Eyis ~ 4.5 MeV, where the observed 
spectrum touches the no-oscillation prediction (see Fig. 9). The best fit 
point lies in the lower region called LMA-I, while LMA-II is allowed only at 
99% CL [16,17]. 

This was followed by the data from the second (salt) phase of SNO [8], 
with better NC detection efficiency. Combining the data from the two phases 
in a global analysis helped to constrain the mass and mixing parameters 
further [18,19]. Fig. 8 shows the results of global fits with phase-1, phase- 
2 and the combined SNO data [18]. The combined fit is seen to allow the 
LMA-II region only at 3a and disallow maximal mixing at 5cr. The most 
important issue at this point was a definitive resolution of the LMA-I and 
II ambiguity. It was shown through a simulation study in [18] that if the 
KL spectrum from 1 kty data continues to favour the LMA-I region, then 
combining this with the global solar neutrino data will rule out LMA-II at 
> 3cr level. 

Recently the KL experiment has published their 766 ty data, whose best 
fit point is indeed in the LMA-I region [20]. Combining this with the global 
solar neutrino data rules out LMA-II at > 3a while sharpening the LMA- 
I region further [21,22]. Fig. 9 compares the KL spectrum with the no- 
oscillation and oscillation best fit predictions. Fig. 10 shows the result of 
combined fit to the KL and global solar neutrino data [21]. The 90% CL 
contour of the global solar fit is shown for comparison. In comparing the 
two 90% CL contours we see that the mass parameter is mainly determined 
by the KL data, while the mixing angle is determined mainly by the solar 
neutrino data. The best fit values along with la errors are 

Am^ = (8.0 ± .6) X 10"^ eV^ sin^ 9 = 0.28 ± .03. (24) 

The best fit value of Am^ corresponds to the 1st oscillation minimum (A = 
{£) = 180km) occuring at Ems — 5 MeV, in agreement with Fig. 9. It is 
evident from eq. (j^^ and Fig. 10 that the solar neutrino oscillation has 
finally entered the arena of precission physics. 

Concluding Remarks: Oscillation analysis of global solar neutrino plus KL 
data has been extended to three neutrino flavours [21]. The resulting mass 
and mixing parameters agree very well with those of the two-flavour analysis 
discussed above. Moreover dropping one of the solar neutrino experiments 
from the global analysis is also found to make little difference to the result. 
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Thus the results are stable and robust. As regards future prospects, the 
precission of Am^ will improve further with accumulation of more KL data. 
One expects some improvement in the precission of the mixing angle from the 
3rd phase of SNO, but not from KL. The main reason for this is the occurence 
of an oscillation minimum in the middle of the v,, spectrum [Ems — 5 MeV). 
This means that the coefficient of sin^ 1Q in the oscillation probability (eq. 
IS} is very small over this energy range. An interesting suggestion made in 
[23] is to reduce the base line length by half (~ A/2), which would mean an 
oscillation maximum in the middle of the z7g spectrum instead. A KL type 
experiment at a reduced base line length of ~ 70 km has been estimated to 
improve the precission of sin^ Q significantly. 

It is my pleasure to thank the organisers of PASCOS'04 for their in- 
vitation and kind hospitality. Le me also take this opportunity to thank 
my teammates Abhijit Bandyopadhyay, Sandhya Choubey and Srubabati 
Goswami. 
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Figure 1: The SSM prediction for the solar neutrino fluxes is shown along 
with the energy ranges of the solar neutrino experiments [2]. 
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Figure 2: Schematic diagram of the effective mass (energy) eigenvalues of 



Ve,ii as functions of the solar electron density. 
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Figure 3: The positions of the MSW triangle, the earth regeneration effect 
and the vacuum oscillation maximum a.t E = 1 MeV are shown along with 
those of the SMA, LMA, LOW and VAC solutions. While the former posi- 
tions scale with E the latter ones are independent of it. 
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Figure 4: The predicted z/e survival probabilities (rates) for the SMA, LMA, 
LOW and VAC solutions. 
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Figure 5: The energy (in) dependence of the SK spectrum [5]. 
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Figure 6: Results of global solar neutrino data fits with (left and middle) 
and without (right pannel) the first SNO neutral current data [12]. 
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Figure 7: The 90,95,99 and 99.73% CL contours of global fit to solar neutrino 
plus the 162 ty KL data. The 99.73% CL {3a) contour of the solar fit is shown 
for comparison [16]. 
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Figure 8: The 90,95,99 and 99.73% CL contours from global fits to solar 
neutrino plus 162 ty KL data with phase-I, phase-II and combined SNO 
data [18]. 
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Figure 9: The 766 ty KL spectrum compared with the no-oscillation and 
oscillation best fit predictions [21]. 
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Figure 10: The 90,95,99 and 99.73% contours from global fit to solar neutrino 
plus 766 ty KL data. The 90% CL contour of the solar fit is shown for 
comparison [21]. 
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